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VALIDATION OF ZERO-ORDER FEEDBACK STRATEGIES FOR MEDIUM-RANGE AIR-TO-AIR 

INTERCEPTION IN A HORIZONTAL PLANE 


Josef Shinar* 

NASA Ames Research Center 


SUMMARY 


A zero-order feedback solution of a variable-speed interception game between two 
aircraft in the horizontal plane, obtained by using the method of forced singular 
perturbations (FSP) , is compared with the exact open-loop solution. The comparison 
indicates that for initial distances of separation larger than 8 turning radii of 
the evader, the accuracy of the feedback approximation is better than 1%. The result 
validates the zero-order FSP approximation for medium-range air-combat analysis. 

This feedback solution is a very attractive candidate for airborne implementation. 


1. INTRODUCTION 


There is little doubt that the differential game approach is the most realistic 
mathematical formulation for air-to-air combat problems. As pointed out in a recent 
review paper by Ardema (ref. 1), however, "...this approach is so complex that results 
to date haVe been disappointing." Some hope for a relative breakthrough in this area 
has arisen with the application of singular-perturbation techniques for nonlinear 
zero-sum pursuit-evasion games (refs. 2-6). This method has created the potential to 
generate approximate feedback strategies for an important class of air-combat prob- 
lems, namely, for medium-range air-to-air interception. (Generally, interception 
ranges are characterized by the firing envelope of the respective weapon systems. In 
this context "medium range" is a loose term covering a domain of weapon delivery 
between 4-20 km, excluding both long-range missiles, such as the AIM-54 Phoenix, and 
guns or other short-range weapons.) 

The approximate solution obtained by the method of singular perturbations seems 
to be a very attractive candidate for airborne implementation. The practical use- 
fulness of the approximation has to be evaluated not only by its feasibility, but 
more importantly by its accuracy. For accuracy assessment of any approximation the 
knowledge of the exact solution is required. In the last decade, several numerical 
algorithms were developed to solve the multidimensional, nonlinear, two-point 
boundary-value problem formulated by the set of necessary conditions of game opti- 
mality (refs. 7-14). Even if it is virtually impossible to prove that the converged 
numerical solution of the necessary conditions is indeed optimal, that is, that it 
satisfies the sufficiency conditions for saddle-point optimality (ref. 15), the only 
available way to assess the accuracy of an approximation is by comparing it with an 
outcome of some iterative computing algorithm. Unfortunately, most of the computer 
programs (refs. 7-12) are inoperative. The presently active computing algorithms 
(refs. 13-14) were also unable to provide the required comparison. Communication 


^Senior NRC Associate on sabbatical leave from Technion, Israel Institute of 
Technology, Haifa, Israel. 
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with the research personnel involved revealed that the algorithms are currently 
limited to the solution of fixed-time problems. Moreover, there is a major problem of 
convergence for the relatively long duration of a medium-range interception (ref. 16). 

Lately, another opportunity has appeared for validating the zero-order feedback 
approximation obtained by using singular-perturbation methods. An open-loop solution 
of the variable-speed aircraft versus aircraft pursuit-evasion game in a horizontal 
plane was published in several recent papers (refs. 17-19). This research effort is 
continuing at Ames Research Center. The open-loop solution method allows a relatively 
detailed comparison with a feedback approximation of the same problem, using an iden- 
tical aero-propulsive model. The comparison has been limited to the already existing 
numerical results. 

The objective of this paper is to report the outcome of the above described com- 
parison and to draw appropriate conclusions. In section 2 the medium-range air-to- 
air interception game in a horizontal plane is formulated and the necessary conditions 
for saddle-point optimality are summarized. In section 3 the zero-order feedback 
solution of the game is developed, using the technique of forced singular perturba- 
tion, summarizing the approach detailed in references 2-5. In section 4 the basic 
concept of the open-loop solution (refs. 17-19) is described, and in section 5 the 
conditions of the comparison and the aerodynamic and propulsion models are detailed. 
Results of the comparison are presented and discussed in section 6 and the conclu- 
sions are drawn in section 7. 


2. FORMULATION OF THE AIR-TO-AIR INTERCEPTION IN THE HORIZONTAL PLANE 


2.1 Equations of Motion 

The geometry of the relative position of two airplanes in a horizontal plane is 
depicted in figure 1. In an interception scenario, one airplane (having a superior 
weapons system) is the pursuer P, and the other is the evader E. The dynamics of the 
relative geometry can be expressed either in polar coordinates. 


R = -Vp cos(Xp -</))+ Vg cos(Xg - 1^); R(to) = Rq 
'P = [-Vp sin(Xp - + Vp sin(Xp - '^)]/R; 'l'(to) = 1'q (2) 


or in a Cartesian coordinate 

system. 




^■''e 

cosxp - Vp cosXp 

9 

X 

rt 

o 

II 
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o 

(3) 

^“''e 

sinXg - Vp sinxp 
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n 

o 

It 

o 

(4) 
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II 

1 
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Figure 1.- Geometry of horizontal interception. 
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Aircraft dynamics for trajectory computations can be adequately represented by a 
point-mass approximation. Inmost studies (refs. 3-6, 17-24), constant weight and 
thrust aligned with the velocity vector are also assumed. This set of assumptions 
leads to 


P ^\w w/p 

= g(Tp - Dp) ; 


(7) 

’e ■ 4 - 1 ), 

- - Dp) , 

w = \ 

(8) 

Xp = {-^) 'S Up = 



(9) 

Xe “ ((^) ^8 Ug ■ 


; Xg(to^ “ Xg 

0 

(10) 


The normalized propulsive thrust and the aerodynamic lift and drag forces can be 
expressed by 

l=n-'l(h.M) V(l) 

5 . 5 . q(h,M) C„(Cl,M)/(|) (13) 

where 5 is the throttle parameter in the range of 

0 < c < 1 (1^) 


where q is the dynamic pressure defined by 


q(h,M) = Y 

and where the relationship between drag and lift is approximated by a parabolic polar 

Cp(Cj^,M) = Cdo(M) + K(M) C2 (16) 


The maximum admissible value of the lift coefficient is limited by aerodynamic 
phenomena (e.g., stall, buffetting, or instability): 


Ct 1 


■■raax 


(17) 


The total lift force is also limited because of structural constraints, 
limitation is usually expressed via the aerodynamic load factor. 


This 


n < n„„„ 
— max 


(18) 
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leading, for horizontal flight, to a bound of the aircraft bank angle: 


|y| < y = cos 
' ' — max 


-1 _L 


( 19 ) 


max 


By combining equations (12), (13), (16), and (17), the total drag-to-weight 
ratio can be expressed by 

D = Dq + n2 (20) 

Dq being the zero lift drag and the Induced drag for level flight (n = 1) 

divided by aircraft weight: 

Do = CDo(M) q/(|) (21) 

D. = K(M) (|)/q (22) 

The nondimensional quantity (^)/q the lift coefficient for level flight, 

(|)/q 4 CLj^(h,M) (23) 


For a horizontal flight, where the speed of sound remains constant, it is convenient 
to replace the velocity as a state variable by the Mach number. Accordingly, the set 
of equations (I)-(IO) can be rewritten in a nondimensional form, using equations 
(18)-(23), as well as a normalized time scale 

T-(f)[t-tol (24) 

and, ( )’, denoting derivations with respect to nondimensional time: 


<“o>P - = 

O 

II 

o 

(25) 


o 

II 

O 

(26) 

X’ = (n| - 1)‘^^/M, = “p 

; Xp(0) = xpg 

(27) 


; Xg(o) = xeq 

(28) 


Here, u , u^ are the normalized controls of the pursuer and the evader respec- 
tively, bounded by 


|up| < 1 ; lugl < 1 (29) 

and the maximum turning rates are determined by the physical constraints (17) and 
(18) . These constraints determine the region of instantaneous maneuverability 
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of the aircraft, as depicted in figure 2. Their intersection determines the so- 
called "corner velocity," or Mach number, 

»c ■ '"maxi)/®-' (30) 


and the corresponding maximum instantaneous turning rate. At these flying condi- 
tions, the turning radius of the aircraft is minimal and is given by 


c — c c 

Or using nondimensional variables, 


r„ 4 V„/x„ = v2/g(n2^^^-l)'/" 




1/2 


(31) 


(32) 


Substituting equation (30) into (32) and using equations (15) and (23) yields 


r 

c 


n 


max 


^"Jax - 


1/2 


Cli (M = 1) 

“Ct (M ) 
^max'^ c 


(33) 


In figure 2 the line of best sustained turning rates, obtained equating maximum 
thrust and drag, is also shown: 


x' = (n2 - l)^/^/M 
s s 

where the sustained load factor is given by 


n 

s 


[(T - Dn)/D.]^/2 

max ^ 1 


(34) 


(35) 


Figure 2 also indicates the admissible region of Mach numbers in horizontal flight. 
The lower bound is a consequence of equation (17) with n = 1, while the maximum Mach 
number can be either "placard** limited (because of structural and aerothermodynamical 
considerations) or determined by the equilibrium in straight flight: 


M 


. < M 

min — 


< M 


max 


^min = fCLi(M = 


1/2 


(36) 


M = min 
max 


M 


limit * 


M ^ 


max 

M 


arg = D. + 


max 




The relative trajectory equations can be also rewritten in nondimensional form 
in polar coordinates as 

R' = -Mp cos(Xp - t|<) + Mp cos(Xg - 'P) ; R(0) = Rq (37) 
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Ip' = [-Mp sin(Xp - i|») + Mg sin(Xg - >|/)]/R ; <i<(0) = (38) 

or in a Cartesian form as 

X' = Mg cosXg - Mp cosxp ; X(0) = (39) 

Y' = Mg sinxg - Mp sinxg ; Y(0) = Y^ (40) 


2.2 Differential Game Formulation 

In this subsection an air-to-air interception engagement is formulated, based on 
the nondimensional mathematical model described previously, as a zero-sum differential 
game. For such a formulation the following elements have to be defined: game space; 

game dynamics; admissible control sets; role determination of the players; informa- 
tion structure; game termination; and cost function. 

For the game dynamics described by equations (37), (38), and by (25)-(28) , the 
six- dimensional game space (R, Mp, M^, Xp, X^) is bounded by 


(41) 


(M . ) < < (M . ) 

min^g — E — min^g 

The angular state variables Xt>> Xt^) unbounded, but periodical; that is, 

r E 

-h 2 tt H l/) 

'I'p + 2 tt E Xp (^2) 

Xj. + 2TT E Xg 

The control set of the players is (C , Up) and (Cp» » respectively, subject 

to the constraints (14) and (29). The roles of the players (pursuer and evader) are 
determined at the outset by the superior weapons system of the interceptor airplane. 
This role determination remains unchanged during the duration of the game. 

The_game starts at some Initial conditions (Rq, i|jq , Mpq, Meq, XPqj XEq) ^ 
moment (t = 0) when both players become aware of each other. For the duration of the 
game, perfect (memoryless) state information is assumed. The game terminates when 
the pursuer succeeds (at the first time) to approach the evader within a distance 
equal to the firing range of its weapon system. This condition is expressed by 

R(tp = d (43) 


R'(t^) < 0 


(44) 
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The terminal values of the other state variables are free. Thus, equation (43) 
describes a five-dimensional terminal manifold, and equation (44) determines the 
"usable part" of this surface. The time of game termination (or in other words the 
"time of capture") is determined, based on equations (43)-(44), as 

tj = min arg {r = 7} (45) 

t Ro>d 

The cost function of the game is this time-of-capture, 

</ = t ^ (46) 

The objective of the purs_uer is to minimize t-; the evader wants (if capture cannot 
be avoided) to maximize t^. 

The solution of this perfect information zero-sum differential game, for each 
set of initial conditions, is a triplet: a pair of optimal strategies p*(*), e*(*) 

and the saddle-point value of the cost function satisfying. 

c/ (P. e*) ^(p*, e*) 4 c/* ^ jT(p*, e) (47) 

This quantitative game (game of degree) has a meaning only if capture, defined by 
equations (43) and (44) , can be guaranteed. An obvious requirement for satisfying 
equation (44) even for a "tail-chase" configuration is 


Mp(tf) > 


(48) 


Analysis of a simplified constant-speed pursuit-evasion game indicates that cap- 
ture can be achieved from ar^r initial condition (ref. 20), even against an evader of 
unlimited maneuverability (rg =_0) , if the ratio of the "capture range," R(tf) =7 
to the pursuer turning radius rp exceeds the value [2(1 - sin l)]^/2. that is, 

d/7p > 0.563 (49) 

In interceptions where "all-aspect" guided missiles are used, this inequality is 
easily satisfied. It can be thus concluded that the medium-range interception engage- 
ment described in the Introduction is a suitable example for the above formulated 
differential game of degr^ee. Such encounters are also characterized by an Initial 
distance of separation R^ much larger than the firing range d. 

The first phase of the solution of such a differential game is to apply the set 
of necessary conditions for an assumed saddle-point optimality, generating a "candl- 
for the solution. Then it has to be verified that the saddle— point inequality 
of equation (47) or some equivalent sufficiency condition (ref. 15) is indeed satis- 
fied. For differential games of complex dynamic structure the verification phase is 
prohibitive; therefore, most investigations (refs. 7-14, 17-19) result only in 
"candidate" solutions. 
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2.3 Necessary Conditions of Optimality 


The necessary conditions of saddle-point optimality for the previously formu- 
lated differential game are expressed by 

min max = 0 (50) 

^P’ “p ^E’ “e 

The Hamiltonian of the normalized game is 


= 1 

+ cos(Xp 

+ 

1 

Mg cos(Xg -'!')] 


+ 

X^[-Mp sin(Xp - 

'(') + Mg 

sin(Xg - 4')]/R 


+ 

^%^^p(^max>p ” 

(Co)p - 

(D.) (1 + ^£(X' 
ip T max p 

u^p)l 

+ 


(®o)e ■ 

(D.) (1 + m2(x' 
ip E max ^ 

u|)l 

+ 


1yt.(x' ) + constraints 

XE max F E 



The multipliers X“, etc., are the gradient components of the cost function ^ 
(assuming that such a gradient does exist) and they have to satisfy the set of adjoint 
differential equations: 


II 

_ 

(52) 

3R 


= 

8j? 

9'!' 

(53) 

II 

r< 

9,3? 
” 9Mp 

(54) 

^Me “ 

■ 9Mg 

(55) 

^Xp = 

9,3? 

9Xp 

(56) 

^XE = 

9J? 
■ 'Xe 

(57) 


The end values of the adjoint variables ^re determined by the transversality condi- 
tions applied at the terminal manifold R(t^) = d. Since all the other state vari- 
ables are free on this surface, 



By satisfying equations (50) and (A4) at t^, we also obtain 


X-(t^) = -1/R (t^) = l/[Mp^ cosCXpj - - Mpf cosCXEf - > 0 (59) 

The optimal control strategies can be obtained, in terms of the state and adjoint 
variables, by applying equation (50). 


* 1 

5p = Y ^Mp] . 

^Mp ^ 0 

(60) 

Y sign XMj;] . 

^ 0 

(61) 

u^ = sat t(u*)p] 


(62) 


where is the unconstrained optimal turning strategy of the pursuer given by 

<<)p ■ (63) 


A similar expression holds for the evader. 


j It can be shown that a singular thrust subarc, requiring XMp = Xw =0 or 
XMp = Xg = 0, is a candidate for this time-optimal game, only if = lM;|^imit^ 
Mg = (Mlimit)g» respectively. Equation (50) also requires that ^ 


2 

3 

3u^ 


* 

Up = Up 


^0 


(64) 


3^J<f 

9u?. 


* < 0 
* 


(65) 


= ^E 


These conditions imply, for unconstrained controls, that Xj^^ < 0 and Xu > 0; 
that is, maximum thrust (5* = 1), for both aircraft. Consequently, if zero^thrust is 
optimal (5* = 0) because of Xj^p > 0 or Xj^^ < 0, it also requires that Up = ±1 
(or Up = ±1) . Conversely, constrained turning rate does not necessarily imply zero 
thrust. 


To generate a candidate solution from some initial conditions of the game, a 
12-dimensional nonlinear two-point boundary-value problem has to be solved. Several 
numerical techniques were developed to solve this problem (refs. 7-14), all requiring 
an excessive computational effort, impractical for a rapid systematical assessment, 
as well as for a real-time implementation. For such purposes an approximate analyti- 
cal solution is preferred. Such an approximation can be obtained by the method of 
forced singular perturbations, as described in the next section. 
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3, ZERO-ORDER FEEDBACK SOLUTION 


In this section an analytical approximation is outlined in order to provide a 
zero-order solution to the original differential game formulated in the preceding 
subsection. The approximation is based on applying the method of singular perturba- 
tions as described in references 2-6. First, a singularly perturbed dynamic model 
has to be defined. Analyzing this model to zero-order leads to a set of low- 
dimensional problems solvable analytically. Based on these simplified solutions, a 
uniformly valid zero-order approximation can be synthesized in a feedback form. 


3.1 Modeling Considerations 

The success of the singular-perturbation approach depends on the time-scale 
separation of the state variables. In this respect, inspection of the set of differ- 
ential equations (37)-(38) and (25)-(28) leads to the following observations: 

1. If the normalized separation distance R is sufficiently large, as assumed 
for a medium-range interception, the turning rate of the line of sight 4^’, given by 
equation (38), will be very slow compared with the turning rates of the airplanes 
given in equations (27) and (28). 

2. Longitudinal accelerations of an airplane (see eqs. (25) and (26)) are much 
smaller (of the order of ±0.5 g or less) than the lateral accelerations (limited to 
6-8 g) used for turning in equations (27) and (28). 

3. Equations (32) and (38) , describing the relative geometry in polar coordi- 
nates, are strongly connected and have to be analyzed on the same time scale. 

Based on Jthese observations the following hierarchy of state variables can be 
established: R , (slowest), , M^ (faster), and Xp , Xg (fastest), leading to a 

singularly perturbed dynamic model. Such a model can be obtained either by proper 
rescaling of the variables, as in reference 2, or by artificial insertion of the per- 
turbation parameter e , as in references 3 and 4. It has been pointed out (ref. 5) 
that the second approach, called also a "forced" singular-perturbation technique 
(FSPT) , leads to a zero-order solution which is identical to the result of scaling 
transformation. Because FSPT seems to be much more convenient for multiple time- 
scale problems, it is the approach used here. 

The FSPT model of the differential game describing a medium-range air-to-air 
interception in nondimens ional variables can be now fomulated: 


R' = -Mp cos(Xp - + Mp. cos(Xg - 4') 


O 

11 

o 

(66) 

t/;' = [-Mp sin(Xp - tp) + sin(Xp - i[^)]/R 


'I'(O) = rl>^ 

(67) 

- <”o)p - <«i)p 11 

u|) 

O 

II 

O 

(68) 

- <“o>E - «1>E ^ 

“ 1 ' 

o 

II 

o 

(69) 
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Xp ■ <Xmax'p “p Xp<0)-Xp„ 

=^'Xe - <x;ax>. “e Xp(0) - Xe. 


(70) 

(71) 


The cost function JT remains unchanged and the Hamiltonian of this singularly 
perturbed dynamic system becomes identical to the original one in equation (51) by 
setting for the adjoint variables of the FSPT model, denoted by a superscript "e," 


^Mp = E^Mp» ^xp “ ^Xp> ^Xg = X 


= *.2 


'XE 


(72) 


This notation leads to the following set of adjoint differential equations, having a 
similar singularly perturbed structure as the state equations (66)-(71): 




3R 


1 r 


^Xmp = - 


EXmp = - 


e^X 


Xp 


e-^Xv 


'h.ye 

^Mp 

9 J? ~ 

^Mg 

9^ 

3Xp 

9j^ 

^Xe 


X^(tf) > 0 


X(f;(t£) ~ 0 

XMp(tj) = 0 


XMg(tf) = 0 


XxgCtf) = 0 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 


When using this reformulation, the other necessary conditions of optimality — 
as equations (60)-(65), which were derived from equation (50) — remain unchanged. 

The zero-order FSPT analysis consist of solving subsequently the following 
lower-dimensional problems: (1) a reduced-order game; (2) an outer boundary layer 

(velocity dynamics); and (3) an inner boundary layer (turning dynamics). 


3.2 The Reduced-Order Game 

Setting € = 0 in equations (60) -(71) and denoting the variables by a super- 


results in the following set of equations: 
(R^)' = -M^ cos(Xp -:{»’')+ cos(Xg - '1^’') 

“^(0) = Rq 

(79) 

= [-M^ sin(xp - sin(xp - ij;^)]/^ 

'I'^(O) = 'I'q 

(80) 
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5p(T ) = (Dq) + (D.) 

r max p ^ p ip 


+ (Di), 


( 81 ) 

(82) 


= 0 


(83) 


u^ = 0 


(84) 


The Hamiltonian of the reduced game is 

= 1 + [-m|^ cos(Xp - ' 1 '^) + Mg cos(Xg - 'i'’^)] 

+ [-M^ sln(Xp - + M^ sinCXg - i|<^)]/^ + constraints (85) 


The constraints include now equations (81)~(84). Setting C = 0 in equations (73)- 
(78) as well, leads to 


(Xi) . - 


djtr 

a"? 


= (X,y/F)(i|;’')' ; x^(t.) > 0 


'J' 




(\l)' = - = xIF(i;;’^)' + (X,y/R^)(R’^)' ; X.y (t.) = 0 


at() 


R 


^ "-f' 


ai^ 




= 0 


= 0 


= -m 5 [xi sin(Xp - i>^) - ^ cos(Xp - 'i'’^) ] =0 
9x; ^ R 


= -M^ [xf. sin(Xg - 4'^) - z| cos(Xp - il'^) ] = 0 

9 Xr. ^ R 


( 86 ) 


(87) 


( 88 ) 


(89) 


(90) 


(91) 


,r „r r r 


The last four equations indicate that the "fast" variables Ml, Xp» Xp* have to 
be considered as additional control variables of the players constrained 5y equatio 
(81)-(84). The last two equations lead to 


y equations 

tg(Xp - = tg(Xg - '^^) = tg 
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implying that 


= (Xp - = (Xp - '1'’^) ± niTT , m = 0, 1 ... (93) 

Substituting (92) into (87) yields 

(Xj)' = (->^ ± M^) [Xl sina^ - (xJ/F) coso’^) = 0 (94) 

and consequently 

“ <”> 

leading to the conclusion that 

= m^TT ; m^ = 0, 1 ... (96) 


The correct values of the free ambiguity parameters m and m^ can be determined 
from the saddle-point optimality of the reduced-order Hamiltonian. Substituting 
equation (95) into (85) and applying 


min 


max 


= 0 


_.r r i- 


(97) 


leads, when combined with equation (67), to the following results 



(98) 


Further application of equation (97), combined with equations (81) and (82), 
yields 


r max p 

K = 

max £ 


(99) 

( 100 ) 


( 101 ) 


with a similar expression for the evader. Moreover, substitution of equations (98)- 
(100) into equations (85) and (97) leads to 


’'i' 


( 102 ) 


In summarizing all these results, it can be stated that the solution of the 
reduced-order game is a tail-chase flown by each airplane at its maximum speed. The 
optimal control strategies are given by equations (98)-(101). These strategies are 
playable (ref. 15) only if that is, if xL > 0, as required by 
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equation (44). Since the reduced-order solution does not satisfy the initial condi- 
tions imposed on the fast variables (Mp, Xp> Xp) boundary- layer solutions are 
needed. 


3.3 Zero-Order Outer Boundary Layer (Velocity Dynamics) 

The first (outer) boundary layer is obtained by performing a time-stretching 
transformation 


= t/e (103) 

on the equations (66)-(78) and setting again e = 0. This operation leads to the 
following equations, where the boundary layer variables are denoted by a superscript 
"o": 


^.0 R».E„ 

dT° " 


dij> ^ „ , o 


dt 


0 

o ^ ^0 


_ _ _ 

^ = ^Xax>p - (Vp - (Vp » H^(0)=Mp„ 


_o 


dx 


o = ^E(^max>E - - (°i>, ’ ^(0) = MEo 


“P 


= 0 




= 0 


dx° 


0 => X^(t°) = aI = 1/(M^ - M^) 


dx: 
-Jl = 


dx 


0 => X°(x°) = X^ = 0 
o ip \p 


dX 


Mp 


dx 


lira X - XM_(t^) - 0 
o Mp f 

T 00 


dX^. 


dx' 


E _ 




lim' ^Me 0 

o 

X -»■ 00 
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(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 
( 112 ) 
(113) 



0 


(llA) 


(115) 


The Hamiltonian of this outer-boundary-layer game is 


.A"” - 1 + cos(Xp - ♦,) + M° cos<x° -♦„)]+ XMp [Cp (T^,)p 

- (P„)p - (Dj)p) + xgp [C° - (i,)p - (Dp)^)I + constraints ( 116 ) 


up 1 p 

Game optimality requires 


yielding 


min max = 0 

pO ° r° ° 

^P’ ^E’ 


o o , 
Xp - Xp 'I'q 


(117) 


(118) 


.p - 2 




Sp ^ 0 

(119) 


(120) 


It is easily seen that this game can be decoupled into two independent optimiza- 
tion problems. It can be also shown that for each airplane there exists an equilib- 
rium point determined by = (Mjnax)p respectively, Mg = (Mmax)r;» coinciding 

with the reduced-order solution. To reach these equilibrium points, which automat- 
ically satisfy the matching conditions, it is required (assuming Mpo < (Wmax)p> 

Meq < (Mmax)^) that 


s'" = 5"" = 1 


(121) 


The values of the adjoint variables and can be determined separately, by 

assuming that the opponent has already reached equilibrium, from equation (117) , in a 
state feedback form: 


iMp = - ^^/f^^max^p " ^°0^P ' ^°i^p^ 

- (^o)e - = V^> 


( 122 ) 


(123) 
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It can be verified that equations (122) and (123) always satisfy equation (117) 
and that they are also compatible with equations (119)-(121). This outer boundary 
layer is therefore an accelerating straight-line "dash" for both airplanes in a tail- 
chase; that is, both velocity vectors are aligned with the line of sight. Since this 
direction does not satisfy the initial conditions Xpq , Xeq an additional boundary 
layer is required. 


3.4 Zero-Order Inner Boundary Layer (Turning Dynamics) 

In this second (inner) boundary layer a new stretching transformation 

= t/e^ = T°/e (124) 

is used and the variables are indexed by a superscript "i." This transformation 
results, setting e = 0, in the following: 



= R° = 

(125) 


■ 1 .o , r , 

^ 

(126) 


o 

II 

O 

II 

o 

II 

(127) 


O 

II 

o 

11 

II 

(128) 

dXp 

dt^ 

tx^(Mpo)]^ax4 ’ xJ(0)=Xpo 

(129) 

dXE _ 

dT^ 

tXp (MEo)lmax “e ’ “ ^Eq 

(130) 


X^=X|=XI=1/[M^-M^] 

(131) 


s' = S° “ S' = ° 

(132) 


- x° 

% ■ ^Mp 

(133) 


X^ - X° 

(134) 

"ip _ 

: llm XxT ' - » 

^Xp 00 

(135) 
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dX 


xe 


dT ■ 


— - ; lira Xxg = ° 


3X 


-)- 00 


(136) 


where the Hamiltonian of this inner boundary layer is 


, 7 ^^ = 1 + xl [-MPq cos(Xp - 'I'g) + Meo cos(Xg - 'I'q)] 


'Mp j4 <Wp - <”0>P - 'Vp II + "Po IC/pPflj 






(x* ) + ^yt 7 (X* ) + constraints 

^max p P aE ^'^max £ E 

The necessary conditions of saddle-point optimality require 


(137) 


min 


max = 0 


^1 1 »i 1 

^P* ^P ^E 


(138) 


leading to 


.i 1 


Sp = 2 [1 - sign = 1 ^ 0 


Mp ■ 


= |[1 - sign xSg] = 1 xSg ^ 0 


(139) 

(140) 


u^ = 


sat[(uh ] 
^ P 


i 

'^E 


= sat[(uj) ] 


(141) 

(142) 


with 


(AK = »xJ/2xSp 4o«>i) 


(143) 


(“?>p = »XJ/2XS, h|„(\)/xU>, 


(1A4) 


From equations (129), (130), (135), (136), (143), and (144) the separability of 
the turning boundary-layer game into two independent optimal control problems can be 
concluded. Based on this Independence, the inner boundary-layer ad joints can be 
eliminated from the Hamiltonian (138) assuming equilibrium of the opponent. Substi- 
tuting equations (143) and (144) separately into equation (138), and using equations 
(122), (123), and (102) results in the following respective feedback expressions 
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(145) 


2 

(xjp) = - Mpo)[l - cosCxp - ° ^^4’ ^^0 ♦ '''o> 

(Xxg)^ = Me^(M^ - Me^)[ 1 - cos(Xg - 'f'o)]/(Mp - = A^(Xg, ^Eq* 


The requirement of a stable solution indicates that the turning has to be oriented to 
reduce the angular differences relative to the line of sight. Consequently it fol- 
lows, observing equations (143) and (144), as well as < 0 and > 0, that 

sign Xxp = -sign Up = sign (Xp - <i<) (147) 

sign X^^ = sign = -sign (x^ - '('g) (148) 


Suiranaitizing equations (141)- (148) and introducing the notion of best-sustained 
turning rate x’> defined in equations (34) and (35), the strategy pair of the turn- 
ing boundary- layer game is 


= sat 


[Xp(Mpg)]g (Mpg[l - cos(Xp - 


^l/2\ 


JXp(MPo)]^ax 


- Mp( 


sign (T|i- - 


Xp) 


= U(x^. Mpq, 


(149) 


u,. 


= sat 


' tXE(MEg)l3 


“ cos(Xp - 




sign (i|;„ - 


Xg) 


= U(Xp. MEq, tl/g) 


(150) 


3.5 Zero-Order Composite Strategies 

Based on the solution of the reduced-order game and the two boundary- layer 
games, a uniformly valid zero-order composite strategy pair can be synthesized. This 
is done for a multiple time-scale differential game (refs. 2-6) by replacing, in the 
control strategies of the fastest (innermost) boundary layer, the initial conditions 
of the slower state variables by their actual current value. (This is clearly a 
feedback solution, where the current state is considered as a new set of initial con- 
ditions.) Such strategies guarantee that as the state approaches the equilibrium 
predicted by the reduced-order solution, the controls also become those required to 
maintain this equilibrium. Moreover, the state feedback form takes implicitly into 
account the eventual (small) variations of the "slow” variables in the boundary layer. 

Comparing equations (139)-(144) with equations (60)-(63) indicates that the 
functional form of the control strategies of the exact optimal and the zero-order 
FSPT composite solutions is identical. The only difference is that in the FSPT solu- 
tion the adjoint variables are approximated by zero-order feedback expressions. 


20 



For the sake of completeness these expressions are summarized here, denoting 
the variables of the zero-order composite solution by a superscript "c": 


= 


= const 


R (m; - M^) 


(151) 


X^ = 0 

'I' 


(152) 


c Mp ■ ^ 1 

XMp = - ^ 7 • -z z — Z 

M_ - (T - D- - D.) 

P E max 0 


(153) 


- M 

c E 




M; - (T - Df, - D.) 

P E max 0 ±' ^ 


-™p - «p 




.1/2 


'Xp - r r 
Mp - 




(D ) [1 - cos(x„ - i^)]^ sign (x„ ~ ’!') 


(154) 


= Xp, <!>) 


2M 2 (m^ - M _ 

M~ tl-cos(x - 

( "e E 


1/2 

l^)]^ sign (Xg - 


(155) 


= Xp. 'P) 

The uniformly valid zero-order composite control strategies are 

«P ^ 




> «P ■ 


(156) 


(157) 


c 

^E “ 


'[(D. + D,)/T ] 

‘ ^ 0 i' max E 


= sat 


[Xp (Mp)] (M^ [1 - cos(Xp - tf')] ( \ 

jsign (ip 


IXp (Mp)]jjjg^ 

U(Xp, Mp. rp) 


. Mp-Mp 


»E ^ 

^ “ ^^limit^E 

S } 


Xp) 


(158) 


(159) 
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(160) 


Ue = 


K (Vis 


[Xg (Mg) ] 


max 


1 = 


[1 - cos(Xg - 'll) 


►^-Me 


1 ) 


sign (i|/ - Xg) 


= U(Xg, Mg. ri>) 


The following qualitative comments can be made about these zero-order strate- 


gies. 


1. Full thrust (5 = 1) is an inherent feature of the FSPT model selected for 
the analysis (if the placard-speed limit is not reached). It is a direct consequence 
of the assumption that speed variations are slower than turning dynamics. In most 
cases this assumption is verified by the solution itself (see section 6.1). However, 
"exact” numerical solutions reveal that for some extreme set of initial conditions, 
such as combination of very high speed and large angular disadvantage (x " ^ ^ + tt) 
zero thrust is optimal. To approximate such situations a different FSPT model, 
described in reference 4, is needed. 

2. The commanded unconstrained (subopt imal) turning rate can be expressed by 

t 

( 

(X')°l .x'(M) (161) 

® ( m’' - M ) 


Relating the commanded turning rate to the best sustained turning performance indi- 
cates the compromise made in the optimization process between fast turning and 
acceleration, both required for successful interception. 

3. The function (1 - cos3)^/^ depicted in figure 3 is of particular interest. 
For small values of 6 it is almost linear. This indicates that for small devia- 
tions from the required heading, the optimal strategy reduces to a proportional con- 
trol with a gain of (x^ [M/(M^ - M)]^/^}. 

4. The speed-dependent gain factor provides an additional Insight. It indi- 
cates that if the current speed is near to its maximum value, that is, if there is 
little need for acceleration, a harder turn can be made. Conversely, if (M^ - M) is 
large, it is better not to lose speed in a sharp turn but to use excess thrust for 
acceleration. 


5. It can be shown that the gain factor is monotonically increasing in the 
vicinity of the maximum speed for level flight, as determined by thrust/drag equilib- 
rium, and that it approaches a finite limit at this speed: 


1/2 

lim Xg [M/CK - M)] = 

M M 


(-1/M)- 




(T 

max 



D 


i 



(162) 


6. If the maximum speed is placard-limited, that is, if 


M 

max 


M 

limit 


< 


M ^ arg 


T 

max 


= D. 



(163) 
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the gain becomes infinite. It means that if this limiting speed is reached before 
the angular deviation becomes zero, the remaining error has to be corrected using 
maximum turning rate. The same phenomenon was observed in the exact numerical 
solution of similar problems (refs. 20, 21). 

7. The feedback expressions of equations (157)-(160) are easily applicable to 
an onboard implementation. They require measurements that are accessible in most 
airplanes (as speed, heading, and line of sight direction). No range measurements 
or estimation of opponent state variables are needed. Such a feedback chart, super- 
posable on figure 2, is depicted in figure 4. 

The zero-order composite strategies of equations (157)-(160) were implemented 
for a digital simulation on the IBM 360/67 at Ames Research Center, using CSMP III 
language . 

Results of this digital simulation, using the zero-order FSPT composite strate- 
gies, were compared with the set of existing exact open-loop solutions described in 
the next section. 


4. DESCRIPTION OF THE OPEN-LOOP OPTIMAL SOLUTION 


The original horizontal Interception game formulated in section 2 is of six 
dimensions. It has been well known that since aircraft dynamics is independent of 
the horizontal flight direction, the order of the system can be reduced to five by 
using relative angular coordinates. In references 17-19 directions are measured 
relative to the line of sight, and the final line of sight direction is used as 
reference. In this particular coordinate system, the necessary conditions of game 
optimality are decoupled into two one-sided optimal control problems, which can be 
solved Independently. This approach enables one to generate extremal trajectories 
for each airplane, independently of its role, of the position of its opponent, and of 
the capture conditions. Such open- loop trajectories are obtained by simultaneous 
retrograde integration of the state and adjoint equations with a given set of end 
conditions. For a game solution, two respective trajectories can be superposed, 
separated by the distance of capture measured along the final line of sight, if the 
conditions for capture are satisfied. The construction of an individual extremal 
trajectory is essentially identical to the method proposed and explored for the solu- 
tion of a time-optimal turn to a point with unspecified final direction (refs. 20, 
21 ). 


In order to obtain the open-loop optimal trajectory for any given initial condi- 
tion, a systematical search has to be carried out in the parameter space of the 
terminal conditions (terminal speed, direction, and time). This effort and its fur- 
ther implementation is reported in reference 22. Because of the excessive numbers of 
computations required to match a given set of initial conditions, it was decided that 
only existing extremal trajectories would be used for comparison. Any point along 
such an extremal can be considered as an initial point of some engagement. 

The open-loop extremals of reference 22 were calculated in a set of "normalized" 
coordinates, different from the one introduced in section 2 and used in the sequel to 
derive the FSPT approximation. Time was not normalized, but distances were divided 
by the speed of sound in order to be consistent with the representation of velocities 
by Mach numbers. Moreover, relative geometry was computed in Cartesian coordinates. 
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For the sake of convenience in making the comparison, the same variables were used 
in the computer program of the digital simulation, implementing the zero-order FSPT 
feedback control strategies (see appendix) . 

The aircraft model used in the simulation (identical to the one used in refer- 
ence 22) and the criteria for selecting examples for comparison are discussed in the 
next section. 


5. AIRCRAFT MODEL AND CONDITIONS OF COMPARISON 


In order to be consistent with the limitations of using only existing results, 
all computations were performed for an altitude of 20,000 ft, implying a = 316.1 m/sec 
and p/po = 0.5326. For the same reason the same aircraft model (an early version of 
the F4-B, used first by Bryson et al. (ref. 23) and later by Bryson and Parson (refs. 
20, 21) served both as pursuer and evader. General aircraft information is summarized 
in table 1, and the aerodynamic and propulsion data at h = 20,000 ft are given in 
table 2 for Mach numbers from 0.6-1. 6. 

TABLE 1,- GENERAL AIRCRAFT INFORMATION 


Weight 


W = 25,000 lb (15,875 kg) 

Wing area 


S = 530 ft2 (49.23 m2) 

Wing loading 


(W/S) = 66 lb/ft2 (322.4 kg/m2) 

Maximum Mach number 

(at 20,000-ft) 


Maximum load factor 


^ax “ 

Maximum usable lift 

coefficient 

Cl = 0.72 

^ax 


TABLE 2.- AERODYNAMIC AND PROPULSION DATA 


Parameter 

Mach number 

0.6 

0.8 

1.0 

1.2 

1.4 

1.6 

So 

0.013 

0.013 

0.014 

0.041 

0.039 

0.036 

K 

0.157 

0.157 

0.180 

0.247 

0.296 

0.343 

"^max/W 

0.491 

0.566 

0.666 

0.780 

0.903 

1.02 


Note: Altitude = 20,000 ft; speed of sound = 316.1 m/sec; 

density ratio = 0.5326. 
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Since the singular-perturbation methods assume a time-scale separation between 
trajectory and airplane dynamics, only open-loop extremals of relatively long dura- 
tion can be considered as candidates for meaningful comparison. As a reasonable 
limit, t^ ^ 50 sec was selected. Respective points on any two extremals can be 
candidates for initial conditions if the capturability condition 

Mp(tp > M^(t^) (164) 

is satisfied. The terminal and initial conditions of the available extremal tra- 
jectories are summarized in table 3. For all cases, the optimal throttle setting is 
maximal (5^ = 1), consistent with the FSPT solution described in section 2. 


TABLE 3.- END CONDITIONS OF EXTREMAL TRAJECTORIES 


Trajectory 

number 

^f 

M 

^f» 

sec 

Xq. 

rad 

rad 

Y ^ 

Xf 

Yf- 

I 

1.5 

1.238 

97.1 

1.330 

1.7 X 10“7 

126.68 

13.21 

II 

1.4 

1.174 

97.1 

3.077 

1.7 X 10“'7 

97.30 

16.48 

III 

1.377 

1.2 

50.0 

9.550 

1.2 X 10"5 

63.78 

4.83 

IV 

1.272 

1.2 

50.0 

2.280 

3 . 5 X 10“5 

47.11 

14.20 

V 

1.222 

0.9 

60.0 

2.917 

5 . 2 X 10"^ 

49.05 

11.31 


^One unit is 316.1 m. 


For the selected initial condition pairs listed in table 4, the best turning 
radii of the airplanes are computed. For any pair of initial conditions, the value 
of (Xq - X^) , computed from the individual trajectories as 

Xq - Xj = Xp(t^) - X^(t^) ( 165 ) 

was kept constant. By varying X^ (the capture range), a large set of differential 
games, with different initial ranges, can be formulated. 

For the constant-speed model of reference 2, it was shown that the small param- 
eter, characterizing the time-scale separation between trajectory and aircraft 
dynamics, is the ratio of the minimum turning radius of the airplane r to the 
initial distance of separation 


£g A r/Rg (166) 

Though in the variable- speed case the technique of "forced" singular perturbations 
(an artificial insertion of a perturbation parameter £ = 1) was used, such a ratio 
still provides a measure of the time-scale separation. The validity of the singular- 
perturbation approach depends on the value of this parameter. For each pair of tra- 
jectories the distance of capture selected in such a way that Eg be in the range 
of 0.1 < Eg < 0.25. For sake of uniformity the turning radius of the pursuer was 
used as reference. 
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TABLE A.- SELECTED INITIAL CONDITIONS FOR COMPARISON 


Example 

number 

Trajectory 

number 

M 

PO 

P0» 

rad 

- a 
r 

P 

^0 

^Eo’ 

rad 


A /V 3 

Xo-Rf 

^f ’ 
sec 

P 

E 

1. 

I 

II 

1.238 

1.33 

8.35 

1.174 

3.077 

7.51 

29.38 

97.1 

2. 

III 

IV 

1.20 

0.55 

7.84 

1.20 

2.280 

7.84 

16.67 

50.0 

3. 

III 

V 

1.20 

0.55 

7.84 

0.822 

0.960 

4.43 

21.23 

50.0 

4. 

II 

V 

1.167 

0.03 

7.42 

0.90 

2.917 

4.41 

28.62 

60.0 


One unit = 316.1 m. 

In the next section selected results are presented in detail and analyzed. 


6. COMPARISON OF RESULTS 


6.1 Qualitative Comparison 

Analysis of the open-loop extremals allows certain qualitative conclusions to 
be drawn. The very nature of the independent extremals exhibits a strong similarity 
with the zero-order FSPT solution. As was pointed out in section 3, the forced- 
singular-perturbation technique leads to a decoupling of the original differential 
game to a ^simple pursuit" game and two independent sets of multiple time-scale 
control problems optimizing an accelerating turn to the instantaneous line of sight. 
The exact open-loop-optimal solution does the same, but relative to the final direc- 
tion of the line of sight. If time-scale separation between airplane turning dynam- 
ics and variations of the relative geometry is a valid assumption, then the current 
direction of the line of sight presents a good approximation to its final direction. 
In such a case, it can be expected that the accuracy of the zero-order feedback 
approximation would be satisfactory. 

The control strategies obtained by the two different methods are vey similar. 
Both require full thrust for time-optimal interception if the maximum speed limit is 
not attained. (In the examples used, the speed limits were never reached.) The 
turning strategies are characterized by a gradually decreasing rate of turn, as the 
flight direction asymptotically reaches the required value. The only difference 
between the exact and the zero-order solution is the reference direction. The larger 
the difference between the instantaneous and final directions of the line of sight, 
the worse the accuracy of the FSPT approximation. 

These qualitative differences are shown by comparing the time histories of the 
control and state variables of the pursuer and the evader for a particular example. 
The initial conditions for this example are depicted in figure 5. In order to obtain 
observable differences, a relative high value of the geometric perturbation param- 
eter, 6g = 0.25, was selected. In figures 6 and 7 the load factor time histories of 
the evader and the pursuer are shown. Figures 8 and 9 depict the time histories of 
the respective flight directions; the line of sight direction is also shown in these 
figures . 
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It can be observed from these figures that the initial load-factor time- 
histories are very similar. The evader, using the FSPT strategy, reaches the actual 
line of sight very nearly at the same time as the open loop optimal trajectory 
closely approaches the final (reference) direction. Its FSPT turning strategy can, 
therefore, be considered as close to optimal. One cannot say the same about the 
pursuer. Since the initial conditions of the pursuer are more favorable than those 
of the evader, the pursuer can reach the line of sight in a shorter time. Because 
the evader has not yet completed its turn, the line of sight continues to rotate. 
Consequently, the pursuer, using FSPT control strategy, has made an unnecessary turn 
(an overshoot) which has to be corrected later. As a result, it can be expected that 
the capture time of the subop timal FSPT game will be longer than the optimal one. 

In figures 10 and 11, Mach number time-histories of the evader and the pursuer 
are shown, respectively; in figures 12 and 13 the initial turning geometry of the 
airplanes is depicted. It can be observed in figure 11 that because of the unnec- 
essary turn, the velocity of the pursuer, using FSPT strategy, is slightly lower than 
the velocity shox^ by the results of the optimal solution, as it has been expected. 
This is in addition to the geometrical disadvantage shown in figure 13. Indeed, for 
this example, the time of capture predicted by the FSPT approximation is about 6.5% 
higher than the optimal value of 97.1 sec. 



1 UNIT = 316.1 m 


Figure 5.- Initial conditions for example No. 1. 
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t, sec 


Figure 7.- Pursuer's load factor time-history. 
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Figure 8.- Evader's turning time-history. 



Figure 9.- Pursuer’s turning time-history. 
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Figure 11.- Pursuer's Mach number time-histoiry . 
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6.2 Quantitative Comparison of Payoff Accuracy 

The purpose of such quantitative comparison is to determine the range of param- 
eters or the domain of initial conditions for which the accuracy of the zero-order 
FSPT feedback solution is satisfactory. If this domain contains a region of practi- 
cal interest, then the approximation can be considered useful. 

The example analyzed in detail in the previous subsection and depicted in 
figures 5-13 was recomputed with different initial and terminal ranges, keeping the 
relative "distance-to-close" constant. In figure 14 the relative payoff accuracy of 
the zero-order FSPT feedback approximation, defined as 




(167) 


is plotted versus the inverse of the geometrical perturbation parameter of the pur- 
suer, defined by 




(168) 


as a ratio of the initial range divided by the pursuer’s best turning radius at its 
initial Mach number, MpQ . 

In figures 15-17, similar results of three additional examples, generated based 
on the available extremal trajectories listed in table 3, are plotted. The initial 
game conditions are summarized for comparison in table 4. These examples are char- 
acterized by the fact that the final Mach number of the extremal trajectory is very 
different from the maximum value. Consequently, the assumption made in the FSPT 
solution that M^ = M is not validated. This discrepancy can result in control 
strategies that devia?i^strongly from the optimal one. In the present examples this 
difficulty was overcome by adjusting the value of M in the respective control 
strategies by a min-max (max-min) parameter optimization technique. The values of 
the adjusted parameters Mp, M^ are noted in the respective figures. All four of 
the last figures (figs. 14-17) show that the accuracy level of 1%, which is highly 
satisfactory for all practical purposes, is satisfied if Eg is of the order of 1/8 
or less. Since in all the available examples the pursuer had less of a turning 
requirement, its suboptimal strategy deviates from the optimal more than that of the 
evader. Consequently, the capture time obtained by the FSPT strategy pair is always 
higher than the optimal. 
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7. CONCLUSIONS 


Results of the comparison, presented and discussed in this report, lead to the 
conclusion that the method of forced singular perturbation provides a valid mathe- 
matical tool for analyzing medium- and longer-range pursuit-evasion engagements. The 
demonstrated accuracy of the FSPT approximation is very encouraging. Taking 1% 
accuracy as a measure of practical usefulness, it follows that the method is appli- 
cable to medium-range air-to-air combat (as defined in the Introduction). 

The turning radius of an airplane at speeds higher than the corner velocity is 
determined by equation (31). This expression gives (approximately), for an example 
of a 6-g-limit load factor and an altitude of 20,000 ft, 

r a 1.5 (km) (169) 

The maximum firing range of advanced air-to-air missiles in a tail-chase situation 
(being of the order of 4-8 km, depending on the altitude) is typically larger than 
the aircraft turning radius. As a consequence, side-stepping maneuvers, which would 
have invalidated the presented FSPT approach, are not optimal. Results of the pres- 
ent comparison show that the required accuracy of 1% is guaranteed by Cg <_ 0.125 
(see figs. 14-17). 

The domain of initial ranges for valid zero-order FSPT analysis can be thus 
determined by 


12 (Mq)| (km) (170) 

Such an inequality is definitely compatible with medium-range interceptions of a real 
operational scenario. 

It should be noted, however, that the validation of the FSPT approximation pre- 
sented in this report is by no means comprehensive. Because of inherent limitations, 
detailed in section 5, only horizontal engagements were analyzed, only a single air- 
plane model was investigated, and the set of initial conditions was limited. More- 
over, the necessity for eventual parameter adjustment requires further investigation. 

Nevertheless, the accuracy of the FSPT method was demonstrated recently in an 
independent study of a closely related problem (ref. 24). In that work, the three- 
dimensional medium-range interception, formulated as a one-sided optimal control 
problem, was analyzed. Such a demonstration could not have been performed for a dif- 
ferential game formulation because of the lack of an operating numerical technique, 
as explained in the Introduction. The currently presented comparison, together with 
the results of reference 24 provides a sufficient basis to encourage further investi- 
gations applying FSPT to air-combat-motivated differential game analysis. 

The attractiveness of the FSPT approach is that it provides a feedback control 
strategy for onboard applications. Because of the simplicity of the implementation, 
it can be incorporated in any future integrated fire and flight control system with 
virtually no additional effort or cost. 
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